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1 Introduction 



Although the dynamics of quarks and gluons at the low energy scale is expected to 
be relativistic and strongly coupled, the simple non-relativistic quark model sucess- 
fully describes the properties of most light mesons (qq) and baryons (qqq). However, 
exception occurs for some scalar particles. Just as stated in Ref. [1], "the features of 
QCD not (all) contained in the (simple) quark model". To describe these overpopu- 
lated scalar particles, noia-qq structures have been assigned to them for about three 
decades since the study in Ref. [2]. They have been regarded as four-quark states [2-5] , 
or molecules composed of conventional particles [6-11] (e.g. KK for / (980) and/or 
a (980)), etc.. Up to now, the puzzle about the nature of these scalar particles still 
remains unsolved. For example, to describe the recent experimental data [12] and the 
more accurate measurement by KLOE Collaboration [13], /o(980) and/or ao(980) were 
regarded as four-quark states [5] or molecular binding of KK [9] and both of them 
lead to results consistent with the experiments. Obviously, further investigation on the 
structure of these scalar particles is necessary. 

On the other hand, more and more overpopulated states (especially those contain- 
ing heavy flavors) have been discovered and confirmed by various experiments [14]. 
Due to the proximity of these particles' masses to those of two lowest lying conven- 
tional particles (carrying certain heavy flavor(s)), one would naturally identify them as 
molecules of conventional particles (see, e.g. Refs. [11,15]). Therefore, it is interesting 
to study whether this picture about these scalar systems is right or not. 

In this paper we will focus on the scalar particles / (980) and a (980). One purpose 
of the present paper is to investigate whether the bound states of the KK system, 
interacting by exchanging various vector particles (p, u>, 0), exist. The other purpose is 
to discuss the extent to which the KK component contributes to the observed particles 
/o(980) and a (980) . 

We choose the Bethe-Salpeter (BS) formalism (in the ladder approximation and the 
instantaneous approximation) as our starting point. The main reason for this is that, 
in comparison with the potential model (used in e.g. Refs. [10,11]), one can include 
some relativistic corrections automatically in the BS equation. 

One may wonder whether the ladder approximation taken for our pseudo-scalar 
system in this paper is suitable. In fact, there have been some works in which the 
legitimacy of the application of the ladder approximation in the BS formalism has been 
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studied, see e.g. Refs. [16-21]. For example, in Ref. [16] it was shown that including 
only ladder graphs in the scalar-scalar system can not lead to the correct one-body limit. 
Furthermore, in the gauge theory, within the ladder approximation gauge invariance 
can not be maintained. To solve these problems, at least crossed-ladder graphs should 
be included [16,17]. More recently, it was shown that the crossed-ladder graphs do 
contribute large corrections to the ladder approximation in some cases [18, 19]. For 
large enough coupling, the contribution from the crossed-ladder graphs becomes even 
more important than that from the ladder ones. 

In our case, the square of the effective dimensionless coupling constant (see Eqs. 
(I44j)-(j43) and (03) in Sect. G2) can be written as g 2 KKV E 2 /(AttM^) which is greater 
than 3q From the naive point of view, for such a large coupling constant, the ladder 
approximation is not legitimate [19]. However, a closer examination shows that there 
is a significant difference between our case and the cases discussed in Refs. [18,19] (see 
also Refs. [20,21]), in which the mass of the exchanged particle is very small. On the 
contrary, the exchanged particles, p, u>, 0, in our case have large masses compared 
with the constituent particles K and K . We will show that the large masses of the 
exchanged particles suppress significantly the contribution of the crossed-ladder graphs 
since factors of the form l/(p 2 — My) from the extra propagators in the crossed-ladder 
graphs lead to extra suppression (in powers of 1/My). Therefore, in our case, the net 
contribution of the crossed-ladder graphs is in fact very small compared with that of 
the ladder graphs (more details are given in Sect. [3]). 

Since the contribution of the crossed-ladder graphs is small in our case the problems 
associated with the ladder approximation, if existing, will not be serious. 

Another approximation we will take is the instantaneous approximation. In this 
approximation, the energy exchanged between the constituent particles of the binding 
system is neglected. This is appropriate if the relativistic effects in the system are small. 
Our calculations (in the ladder approximation and the instantaneous approximation) 
show that both the iso-scalar and iso-vector bound states of the KK system with small 
binding energy exist. This shows that the binding of the constituent particles is weak, 
hence the exchange of energy between them can be neglected. 

However, regarding these bound states as the observed particles, / (980) and a (980), 



* In our case, the binding energy is of order 0(1Q 1 ) MeV, the total energy of the binding system 
is E w 2Mk , and the coupling gxxv is about 3 when V = p,uj and about — 3\/2 when V = <fi. 
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one can find that the decay widths of these bound states are too small to explain the 
experimental data. In other words, while the bound states of KK do contribute to the 
observed scalar particles, they themselves can not describe the full properties of these 
particles. 

The remainder of this paper is organized as follows. In Sect. EJ we review the BS 
formalism for the system of two pseudo-scalar particles and discuss the normalization 
condition of the BS wave function. In Sect. EJ we discuss the bound state equations 
for the KK system in detail. The decays of the KK bound state to tttc and ttt7 final 
states are discussed in Sect. HL The numerical results are presented in Sect. EL The 
final section is reserved for some discussions and our conclusions. 



2 The Bethe-Salpeter formalism 

In this section we will review the general formalism of the BS equation and derive 
the BS equation for the system of two pseudo-scalar particles. We will also derive the 
normalization condition for the BS wave function. Let us start by defining the BS wave 
function for the bound state \P) of two pseudo-scalar particles as the following: 

X P ( X ^ X 2) = (0\T(f> l (xi)<f>2(x2)\P} = e- iPX x p (x), (1) 

where (f>i(xi) and (j> 2 (x 2 ) are the field operators of two pseudo-scalar particles, respec- 
tively, P denotes the total momentum of the bound state, and the relative coordinate 
x and the center of mass coordinate X are defined by 

X = 771X1 + 7/2X2 , x = %i - x 2 , (2) 

or inversely, 

xi=X + 7] 2 x, x 2 = X-r]ix, (3) 

where rji = vrtij [vn\ + 777.2) ,mi(i = 1, 2) is the mass of the i-ih constituent particle. The 
equation for the BS wave function can be derived from a four-point Green function, 

S(x 1 ,x 2 ;y 2 ,y 1 ) = (0|T 1 (x 1 )0 2 (x 2 )(0i(z/i)02(y2)) t |O) . (4) 

To obtain the BS equation, we express the above four-point Green function in terms 
of the four-point truncated irreducible kernel K , 

S(xi, x 2 ; y 2 , 7/1) = £(0) (xi, x 2 ; y 2 , yi) 

+ J d 4 u 1 d 4 u 2 d 4 v 1 d 4 v 2 S (0) (x 1 ,x 2l u 2 ,Ui)K(u 1 ,u 2l v 2 ,v 1 )S(v 1 ,v 2l y 2 ,y 1 ) , (5) 
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where S(p) is related to the forward scattering disconnected four-point amplitude, 

S {0) (x 1 ,x 2 ;y 2 ,yi) = A 1 (x 1 ,y 1 )A 2 (x 2 ,y 2 ) , (6) 
where A,(xi, yi) is the complete propagator of the i-th particle, 

A t (x,y) = <0|T^(^)^(2/)t|0> = J e -*<«-i0 A i (k,m i ) . (7) 

From Eqs. P]) and ([5]) one can derive the following BS equation for the bound state of 
two pseudo-scalar particles: 

X P (x!,x 2 ) = J d A u x d A u 2 d A v x d% 2 5(0) (xi, x 2 ;u 2 ,u 1 )K(u 1 ,u 2 ;v 2 ,v 1 ) x p {vi,v 2 ) , (8) 

or, by inverting S(p), 

J d i yid i y 2 S^ ) {xi,x 2 ;y 2 ,y l )x p {yi,y2) = J d 4 Vxd% K(x 1 ,x 2 ;v 2 ,v 1 ) x p {vi,v 2 ) . (9) 

In this paper, we will investigate the BS equation in momentum space, in which the 
BS wave function is obtained as (using Eq. (CQ)) 

XM>P*) = J d'x x d A x 2 e^^xM,x 2 ) = (2n) A 5( Pl +p 2 - P) Xp (p) , (10) 

where p = r] 2 pi — r]ip 2 is the relative momentum and Xp (p) = J d A x e ipx x {%) ■ The 
Fourier transformation of the four-point Green function in Eq. (j3J) reads 

S( Xl ,x 2 ;y 2 , yi ) = J ^^f^ e -^-^ S(p,p', P, P>) , (11) 

with S(p,p',P,P') = (27r) 4 5 4 (P — P')Sp(p,p'). Similarly, for the irreducible kernel we 
have 

K(xu x 2 - y 2 , Vl ) = J ^^t^' e-vx^-Wy K(p, p', P, P>) , (12) 

with K(p,p' , P, P') = (27r) 4 5 4 (P — P')Kp(p,p'). The relative momenta and the total 
momentum of the bound state in the equations are defined by 

P = V2P1 ~ V1P2 , p' = V2P[ ~ ViP 2 , p = Pi + Pi = Pi + P 2 , (13) 
or inversely, 

p 1 = rj 1 P + p, p 2 = r] 2 P-p, p' 1 =r] 1 P + p', p' 2 = j] 2 P - p' . (14) 
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We must note that the constituents of the bound state can not be on-shell, otherwise 
the bound state is not a really bound state. Consequently, p\ ^ mf (and similar for 

Then, the inhomogeneous equation (J3]) in momentum space reads 

/ (0[ / p(P^)+^(P^)]^(^P / ) = (2 7 r) 4 5( J 9- J5 / ), (15) 

where Ip(p,k) = — (27r) 4 <5 4 (p — k)A^ 1 (pi, mi)Aj 1 (p2, ^2) • The BS equation ([9j) for 
the bound state in momentum space takes the following form: 

[lp(p,k)+-K P (p,k)] Xp {k) = 0. (16) 

This is a homogeneous equation for the BS wave function. 

From the BS bound state equation, Eq. (Q in coordinate space or Eq. ffl6|) in mo- 
mentum space, we can see that the BS wave function satifies a homogeneous equation. 
Therefore, its normalization can not be determined from the bound state equation. To 
obtain the correct normalization of the BS wave function, following Ref. [22], we start 
by considering the contribution of the bound state with P° = Ep to the four-point 
function S. Let us first isolate the contributions from some possible bound states. 
Consider the case with Hiinja^, ^2} > maxjyj, y®} and insert a complete set of states 
into the four-point Green function, we have 



S(x 1 ,x 2 ;y 2 ,y 1 ) = J](O|T0 1 ( 2;i )0 2 (x 2 )|P}(P|T0 2 ( 2/2 )t0 1 ( 2/l ) t |O) 



p 



min{a;J,a;§}>max{yJ,j/§} 



(2tt) 3 

Furthermore, the requirement minla^, x®} > max-f^, y®} can be described by a theta- 
function, 

e(x°-Y° + V -^(x°-y°)-^-^ . (18) 

Using this representation and the contour-integral definition of the theta-function, Eq. 
(|17p can be written as 

= i J g ^^^xMxM po _ l Ep + - 

x e -i(P°- El ,)l( V2 - m )(x»~y°)-\x°\-\y°\}/2 _ 
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Hence, near the pole at P = Ep, we have 

Sp{p,p') = ^ p — — X P (p)Xp(p') + terms regular at P° = E P . (20) 

P u — Ep + ze 

Now, define an auxiliary quantity: 

QpM = j (0i( P °-^)^ fc )^o [Ip(*,i/)+^p(*,pO] • (21) 

For convenience, we can imagine the arguments of Kp, Ip, Sp , and Qp as matrix 
indices and write the above quantity in a compact form: 

Q P = (P° - Ep)S P -^ [I P + K P ] . (22) 

In terms of this notation we can also rewrite Eqs. (1151) . (1TB1) . and (1201 as 



Sp[/ P + Kp] = 1, (23) 
[lp + Kp} Xp = 0, (P° = E P ), (24) 



lim (P u -P P )Sp = i Xp X P - (25) 

Using the above equations and operating Qp upon \ P we have the normalization 
condition for the BS wave function, 

1 j d wf^ {p) JP^ [lr(P,P')+Kp(P,P')] X P (P) = 1 , P° = Ep. (26) 

The BS equation ( fl6l) is very complex. Without approximation we can not even 
write down the irreducible kernel and the propagators of particles explicitly. Since the 
binding of the KK system is weak we use the so-called instantaneous approximation: 
Kp(p,p') = Kp(p, p'). Furthermore, the propagator is set to have the form of the free 
one. Then, the BS equation (fl6l) becomes 

- (P? " m\){ P l - ml) Xp ( P ) = j ^K P (p, p') Xp (p') . (27) 

Now, we divide Eq. ( 1271) by the two propagators on both sides and then perform the 
integration over p° and p'°. Then we have 

E 2 - (E 1 + E 2 ) 2 ^ if d 3 p' 



(Ei + E<i)l E\E\ 



X» =\J ^s^p(p,pOXp(pO, (28) 



where E { = \/p 2 + m 2 , E = P° , and the equal-time wave function is defined as 

X» = j dp° Xp (p). (29) 
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While deriving Eq. (1281) we have used the following result (in the rest frame of the 
bound state): 



dp 



—m 



(30) 



(pi - m? + ie) (pi - ml + ie) — (E 1 + £ 2 ) 2 

which can be obtained by choosing a proper contour. For convenience we define the 
following potential: 

i 



V(p,p') 



E\E2(E\ + E2 
Then, the BS bound state equation can be written as 



K P (p,p'). 



(31) 



E 2 



X P {P) 



tzy 

(2^)3 



v(p,p')Xp(p') 



(32) 



L(£i + £ 2 ) 2 

Eq. (I3"2l will be the starting point in our later numerical calculations. 

For later convenience we also write out \ p (p) in terms of x P (p) • From Eqs. (1271) 
and (1281 we have 



x p (p°,p) 



(pi —m\ + ie)(pl — m 2 , + ie) J (27r) 4 
i 1 



' /V K P (p,p')x p (p') 



E 2 — (Ei + E2) 2 
7T (p 2 - m\ + ze)(p| - m 2 + ie) (£1 + E 2 )/E l E 2 

where p\ — m\ = (r/iE + p ) 2 — E\ , p\ — m| = (772-E — p ) 2 — -E" 2 . 



XplPJ 



(33) 



3 The bound state(s) of the KK system 

In this section, we will study the possible bound state of the KK system. The lowest 
lying particles with the strangeness numbers ±1 form two isospin doublets: (K + , K°) T 
and (-K, K~) T , where the superscript T denotes transpose 1. One can gather them 
into two fields, K\ and K2 , which have the following expansion in momentum space 

Kl = Sik^E; [a ^"" + < - e '" ) ' (34) 
R ' 2 = / W?^M {aK ° ei " + a *° e *° ] ' (35) 



t The conventions about the isospin multiplets used here and in the following are the same as those 
used in e.g. Ref. [23] . 
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where E p = a/p 2 + m 2 K is the energy of the particles and we omit the effects of isospin 
violation so that the masses of the kaons are the same. These two fields can be grouped 
into an isospin doublet K = (Ki, K 2 ) T , which furnishes the fundamental representation 
of the isospin group SU (2) f . 

The KK system has isospin 1 or 0. The iso-scalar bound state can be written as 

\P) = -L\K + K- +K°K°) , (36) 
v 2 

and the three components of the iso- vector states are 
|P} 1)0 = -L \K + K~ - K°K°) , \P) 1>+1 = - \K + K°) , |P)i,_i = \K-K°) . (37) 

Let us now project the bound states on the field operators K\ and K 2 . From Eqs. (1341) 
and we have 

(0|T{^(x 1 )^(a; 2 )t}|P) /i/3 = q^ 3) x p (7) (^i^2), (38) 

where x P 1S the common BS wave function for the bound state with isospin / which 
depends only on the state \P)i,i 3 (as will be shown later, the BS wave function only 
depends on / but not on 7 3 ) but not on the concrete field contents. The isospin 
coefficients C l t l T s for the iso-scalar state are 

CSo)=Cgo) = else =0, (39) 

and for the iso- vector state we have 

ClU) = -Cg 0) = V V2 , C7g +1) = -1 , Cg_u = 1 , else = . (40) 

Now consider the kernel. The BS equation ( |271) for the bound state can be written 
down schematically, 

Ar'A-^x/) = Kf'c^x^ , (41) 

where Ai^ are the propagators of the constituent particles. Then, from Eq. (J39l) . for 
the iso-scalar case, we have (take ij = 11 as an example) 

Ar x A 2 V M^V^p' 2 V 0) • ( 42 ) 

Similarly, for the iso-vector case, taking the ^3 = component as an example, we have 

A^A- V 1,0) = (^p 1 ' 11 - ^P* 22 )X P (1 ' 0) ■ (43) 
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From the above equations, we can see that if the iso-scalar bound state exists one can 
not ensure the existence of the iso-vector bound state, and vice versa. 

The interactions among the kaons and vector particles, p, uj, <p, at the level of 
hadrons are described by the SU (3) v x SU (3) A chiral dynamics. The relevant interac- 
tion vertices are (see e.g. Ref. [11]) 

Ckk p = ig KK K\T-p»)d,K + C .c, (44) 
Ckku = is^/tM^ + cc, (45) 
C KK4 > = ig KK K\d,KW + C .c, (46) 

where c.c. is the complex conjugate of the first term and g (V can be p, lu, and </>) 
are the coupling constants which can be related to g p7T7T in the SU(3)f limit, 



9 KKp = 9^/2 , g KKuj = g^/2 , g^ = -g p7T7T /V2. (47) 

The p7m coupling is determined by g p7r7r = M p j (y/2f w ) ~ 6 [24] , where M p is the mass 
of p and f n is the decay constant of the pion. 

From the above observations, at the tree level, in t-channel we have the following 
kernel for the BS equation in the so-called ladder approximation 0: 

K(Pi,P2;P2,p'i,M v ) = -z(2tt) 4 5 4 (p / 1 +p 2 -px -p 2 ) 



.2 (pi + p'i) ■ (P2 + p' 2 ) + ( pI-p?) (pI - p?)/m$ 

\2 _ AJ2 



X ° l9 ^ ( pi _ p /)2_ M 2 > ^ 



where c T is the isospin coefficient: c = 3,1,1 and c x = —1, 1, 1 for p, u, , respectively. 
These results are consistent with those in Ref. [7]. In Eq. fj4*8l) we have used the 
following propagator for a massive vector meson: 

A^(p, M v ) = 2 ~\ A 2 (g»u - P^Pu/M^) . (49) 

P My 

From the above analysis, we can see that the BS wave function depends only on the 
isospin / but not on its component ^3 . This is because we have only considered strong 
interactions which preserve the isospin symmetry. Therefore, we will omit the -Z3 label 
and write \ ^ simply as \ from now on. 



-yy 22 11 11 12 12 11 11 

* When the exchanging meson is p , we have K P ' (p) — 2K P ' (p) and K P : (p) = —K P ' (p) 

11 22 12 12 11 11 

when the exchanging meson is u>, we have K P ' (ui) — and K P ' (uj) = K P ' (to) ; the case for (j. 

is the same as that for uj . 



In the following discussion and calculation in this paper, we will take the kernel (148p 
as our starting point (this is the so-called ladder approximation) and will not consider 
those non-ladder (e.g. crossed) graphs. Before proceeding, let us discuss briefly the 
contribution of the crossed-ladder graphs. Our aim is to estimate the ratio of their 
contribution to that of the ladder graph. To make the calculation of the (4-th order) 
crossed graph tractable, the following simplification are taken: terms carrying p 2 , p' 2 , 
p ■ p', etc. are omitted in the calculation. Since we will stay in the rest frame of the 
binding system, we also set p ■ P and p f ■ P to zero (the instantaneous approximation). 
These approximations will be appropriate if the half width of the BS wave function is 
small enough compared with the masses My and Mk . From the numerical results (in 
Sect. E]), we can see that the half width of the BS wave function is about 0.1 GeV, which 
is indeed very small compared with the masses My and Mk . The physical picture for 
this approximation is that the configurations with small momenta are dominant in the 
model. 

Taking these approximations, we calculate explicitly the contribution from the (4- 
th order) crossed graph. We take the case where the exchanged particle is p as an 
example, other two cases with uj and <fi as exchanged particles give similar results. We 
work with the dimensional regularization method and the minimal subtraction scheme 
while calculating the crossed graph. The numerical results show that the ratio of the 
contribution from the crossed graph to that from the ladder (2nd order) graph has the 
following form: 

u 2 

0.12 -0.01 In p 9 , (50) 
lGeV 2 V ; 

where \x is the renormalization scale. Since only one (but not all) higher order graph 
is calculated, the result depends on this renormalization scale. It is natural to take \x 
to be around 1 GeV, which is the scale of chiral symmetry breaking. From this result, 
we can see that the ratio of the contribution from the crossed graph to that from the 
ladder one is less than 15% (in the case where ui is the exchanged particle the result is 
almost the same while in the case where <ft is the exchanged particle, we have the ratio 
« 25%). 

For comparison, let us reduce manually the masses of the exchanged particles to 
0.15Mk, which is the case discussed in Refs. [18,19]. We find that this ratio will rise 

2 

to 300% or even more: 3.84 + 3.38 In 1( ^ v a • This result is consistent with those shown 
in Refs [18,19]. 
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Now, let us consider only the ladder approximation and proceed by taking the 
instantaneous approximation, p° = and p'° = in the kernel (118]) . and stay in the 
center-of-mass frame of the bound state, P = . Then, the potential in Eq. ([3"T]) due 
to the exchange of a vector meson V becomes (using Eq. (j!4p ) 



V^(p,p';M v ) = Cj £/(p,p';M 



v 



-g 2 KKV (p + p') 2 + 4f]if] 2 E 2 + (p 2 - p' 2 ) 2 /M 2 
E 1 E 2 (E 1 + E 2 ) (p-p') 2 + M 2 



In order to describe the phenomena in the real world, we should include a form factor 
at each interacting vertex of hadrons to include the finite-size effects of these hadrons. 
For the meson (qq) case, the form factor is assumed to take the following form [7]: 

F(k) = !kn^ k=p - p '- (52) 

where A is a cutoff parameter which will be adjusted to give the solution of the BS 
equation. At the lowest order, the BS equation includes F 2 in its kernel, i.e. V — > V-F 2 . 
The most important term in the numerator of Eq. (jBTj) is Ar\\r\ 2 E 2 . Other terms are 
small since the momenta of the constituent particles of the binding system are small. 
After transforming into the form in coordinate space, similar to the case in Ref. [11], 
one can see that the potential is in fact a Yukawa-like potential (the sum of a Yukawa 
potential and several derivatives of the Yukawa potential). 

Then, for the bound state of the KK system, the BS equation ( |32l) becomes 

1 f d 3 p' 



E 2 

- 1 



x/HIpI) = \ J ^V«\p,p')F(k)%W(\p'\), (53) 



{E 1 + E 2 y 

where the effective potential is (depending on isospin I) 

KSW') = 3tf(p,p^M,) + l7(p,p^M w ) + E/(p,p^M,), (54) 
= -t7(p,p , ;M p )+t7(p,p';M &) )+E7(p,p , ;M^). (55) 

If we are interested in the ground state of the BS equation, the corresponding BS wave 
function is in fact rotational invariant, i.e. X (p) depends only on the norm of the three 
momentum, |p| . Therefore, after completing the azimuthal integration, the above BS 
equation becomes a one-dimensional-integral equation, which reads 

X P (/) (|P|) = J d\p>\ V«(|p|, \p'\)x p V\\p>\) , (56) 
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where (|p|, |p'|) are one-dimensional effective potentials: (|p|, |p'|) = 3L r id(M p ) + 
U ld (M w ) + C/ ld (M^) , and V^flpl, |p'|) = - ^ ld (Af p ) + U ld (M„) + U lA {M^) with 

E 1 + E 2 



u ^ M v>- A(2n) 2 El E 2 [E 2 -(E 1 + E 2 ) 2 ] |p 

where 

\4 = -4|p||p'| (2A 2 - M 



3 ) 



2A 2 + 2(2r ]l r ]2 E 2 + |p| 2 + 


Pf) + (|P 


2 _ 


IP' 


2 ) 2 /M 2 


[2A 2 + (| P 


+ | 


P'l) 2 ] [2A 2 + (|P 


- Il 


P'l) 


2 ] 



^2 = 

v 3 = - 



| P | 2 - |p'| 2 ) 2 /M 2 



v 



In 



M v + (|p| + |p'|) 2 



My + 2(2r] 1 r] 2 E 2 + |p| 2 + |p' 
M 2 + 2(2^ 2 £ 2 + |P| 2 + |pf) + (|P| 2 - |P'| 2 ) 2 /M 2 ] In ^ ± M ± jg|l . 



M 2 + (|p|-|p'|) 2 ' 
, ^2A 2 + (|p| + |p'|) 2 



4 The decay width of the i^Ti^ system 

To find out the bound states of the KK system, one only needs to solve the homo- 
geneous BS equation. However, when we want to calculate physical quantities such 
as the decay width we have to face the problem of the normalization of the BS wave 
function. In the following we will discuss the normalization of the BS wave function 

x p (IpD- 

Substituting the relation between x P (p) an d X p (|p|)j Eq. (1331) . and Eqs. (13T1) 
(I5ip into the normalization equation (1261) one arrives at the following normalization 
equation for \ (|p|) (after carrying out some po-i n t e g ra ti° ns with proper contours): 

-^/I^IPI)] 2 *-^/ ^?,'"(|P|)^'"(|P'I)F^« = 1,(57) 



where F( 



2A 2 -M 2 



2A 2 +(p-p') 

= -H{M p ) + H{M U ) + H(M 4> ) with H(M Y 



is the form factor, = 3H(M p ) + H{M U ) + H{M+) and 



and 



R 



-E 



(p-p') 2 +M4 

- 2E\E\ - El){E xm - E 2V2 ) + E 4 (E lVl + E 2 r] 2 ) 

+ (Ef - E 2 2 ){E 3 lVl + 3£i-E£t7i - 3E 2 E 2V2 - E^) 
2E 1 E 2 [E 4 + (El - E 2 2 ) 2 - 2E 2 (El + E 2 )] ' 
After completing the azimuthal integration in Eq. (|57p we have 

-^/%IIpI 2 x p (/) (IpI) 2 ^ 
- ^ // d \p\ d \p'\ IpIIp'I x p {/) (IpI)x p (/) (Ip1) r« = i 
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(5* 



where T (0) = 3T(M P ) + T{M UJ ) + T(M+) and T« = -T(M P ) + T(M W ) + T{M+) with 

r(M v ) = 4, 



2A- - M 2 2A 2 - My 



2 



2A 2 + (|p| + |p'|) 2 2A 2 + (|p| - |p'|) 

2A 2 + (|p|-|p-|) 2 M* + (|p|-|p1) 21 
2A 2 + (|p| + |p'|) 2 M 2 + (|p| + |p'|) 



for each vector meson with mass My and coupling g KKV . 

If the wave function obtained in the previous section (which will be calculated 
numerically in the following section) dose not satisfy this normalization equation but 
gives some constant c 2 7^ 1 for the expression on the left hand side of Eq. (1551) . one need 
only make the replacement \ (|p|) ~~ > X p (IpI)/I c I to ensure the correct normalization 
of the BS wave functions. 

If the molecular binding is dominant in the KK system, then the possible bound 
states of the KK system are most likely related to the two particles which are denoted 
by / (980) and a (980) in the review of PDG [26], since they are just below the thresh- 
old of the free KK system and, up to now, can not be assigned with the common qq 
structure. One possibility is that they are (mainly) molecular states of other conven- 
tional particles, e.g. K and K. However, to identify the possible molecules of KK with 
these scalar particles, we should also identify other properties (other than the binding 
energy) of the molecules with those of the scalar particles measured by experiments. 
Among these properties, an important one is the decay width of the bound state. 

Now, we will proceed to study the decay widths of the KK bound states and 
compare them with those of /o(980) and ao(980) and see whether the assignment of 
the molecular states with them is suitable. Since the dominant decay channels of 
/o(980) and ao(980) are tttt and tjtt , respectively, we will study the decay widths of 
the above bound states into tttt and t]tt. The relevant interaction vertices are (see e.g. 
Ref. [11]) 



C"kkk* — ig 



■kKK* 



duK*(T ■ tt)K*" -tf {r-d^K^ +c.c. , (59) 



CnKK* = ig vKK * [d^K^-n - tfK*»dtf\ + cc. , (60) 

where cc. denotes the complex conjugate of the previous terms. The coupling con- 
stants are related to g pn7T in the following way: 

9 wKK * = gpnir/2, g vKK * = -VSgpTm/2. (61) 
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The differential decay width of the bound state can be written as [26] 

dr = ^ M ^ dn ^ (62) 

where |q| is the norm of the three- momentum of the particles in the final state in the 
rest frame of the bound state. M. is the Lorentz invariant decay amplitude of the 
process. The lowest order decay amplitude can be written as (the decay to ni] will be 
considered later on) 

(7r a ( gi )7r 6 (g 2 )|^ J d A xd 4 yT {C 7rKK 4x)C wKK ,(y)}\P) 



y/2E 7T a2E 7T t J (2tt 



\4 



k=q+p+{rii-ri2)P 



k=g—p 



x (2tt) H\P -qt- q 2 ) , (63) 

where (i = 1, 2) is the momentum of the z-th particle in the final state and E n a = 
i/q 2 + m\ a , E n b = a/iI + m ^i> ■ The coefficients in Eq. ( l63l) are the isospin 
factors which have been given in Eqs. (1391) and (HOl) . g = 772^1 — 771^2 which is not the 
relative momentum of particles in the final state (note that 771 and 772 are defined as 
rji = vrtij (mi + m 2 ) , and mi and m 2 are the masses of the component particles of the 
bound states but not of the final states). In deriving the above equation the following 
propagator for the vector kaons has been used: 

(0\T{Kr(x)K*» (y)t}|0> = J 0^e-^A, u (k,M K *)5 tJ , (64) 

where i and j are isospin indices. 

The Lorentz-invariant decay amplitude of the KK bound state to tttt is then 

K) ] = ^C,^ J ^{±F(MfA^M K ,Y\k + 2q 2 Y(k-2qi)\ = ^ 

+(k^q + p + (vi-m)P)}x P [I \p), (65) 



where "+" and "— " in "±" are for the iso-scalar and iso-vector channels, respectively. 



For the iso-scalar channel, from Eq. 



Since the iso-scalar tttc final states reads 







we have (r a T%C 3 l = -j=Tr(r a r b ) = y/2& 



ab 



I TTTT) (o,0) = 7= |vr + 7T + 7T 7T+ + 7T 7T ) = -= | 7rV + rV + 7T 3 7r 3 ) , (66) 

v3 ' v3 



§ For the states in isospin multiplets, we use the same conventions as those in Ref. [23] . That is 
to say, we have \ir, ±1) = T-K*) = T'^fl 71 ' 1 ± i^ 2 ) and \n, 0) = | tt 3 > . 
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we have cZ7^ = —y/6. For the iso- vector channel, using Eq. (HOI we have cfe = 2 
since the iso-vector ittt final state is 

l^\ l^±^0 ,^0 +\ ^ 1^0 1 10 I •_0_2 •_2 0\ 

1^/(1^1) = /= P ^ — TTTT ^ = — 1 7T 7T — 7T7T ± 27T 71 =|= 271 71 J , (01) 

. , X I _| _|_ \ % I -1 Q Q -1 \ , , 

I 71 " 71 ") (1,0) = /= P" 7T _ — 71~71 ) = — r= 71 71 — It It ) . (68) 

v2 ' v2 

Note that this isospin coefficient is independent of the component J 3 (notice that our 

convention for the kaon state is different from that in Refs. [7, 10]). This iso-vector 
final state is anti- symmetric, so there is in fact no S-wave 7171 final state with 1=1. 

Now, let us turn to iri] final state. The lowest order matrix element for the decay 
of the KK system into 7177 is 

(7T a (q 1 )r ] (q 2 )\2^ J tfxtfyT {C %KK *(x)£ vKK *(y)}\P) 

= w£^ (r% ^ r ^{ Fm f^ v{K M K ,) % \k + 2q 2 nk~2 qi r 

yf 2E 1T a2E rl y 1 J < k=q-p 

+ (k^q + p+( Vl - V2)P)}x p [I \p) (2tt)V(P -q x - q 2 ) . 
Then, the Lorentz-invariant decay amplitude is (only the iso-vector channel contributes) 
M (i) ~ tc (i) 9„ KK »9 vKK *v2L 

X / -^AF{M) 2 ^Ak,M K *)i 2 {k + 2q 2 y{k-2q 1 y 

J (27TJ 4 I k=q-p 

+ (k^q + p+( Vl - V2)P)}x p (1 \p) • (69) 
The components of nr} are 

N)(i,±i) = TK^} = T-j= I (tt 1 ± in 2 )r)) , kv)(i,o) = I A) = k 3 v) • (70) 

Then, from Eq. (HOI , we have cffl = \[2 which is again independent of 1% . 

In the calculation we stay in the rest frame of the bound state and hence P = (E, 0) . 
In this frame the momenta of the two particles in the final state can be taken as: 
qi = (E[, q), q 2 = (E' 2 , -q) . Therefore, q = r} 2 q\ - r\\q 2 = (r] 2 E[ - rjiE' 2 , q) . When the 
final state is nir , E[ = E n a and E' 2 = E w b while when the final state is irr] , E[ = E n a 
and E 2 = E rj . To calculate the amplitude, we first carry out the azimuthal integration 
of the spatial part of p, the result having the following structure: 

/ nW^fr, M K *)i\k + 2q 2 r(k - 2 qi y X /\p) 

/oo poo 
dp° d\p\\p\ 2 f(p ) Xp {I) (±p ,\p\), (71) 
■oo JO 



(27T) 
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where 

ftf) = (Mk - 2A2 ) 2 

J ^ ; (-2p°q° + |p| 2 + |q| 2 + s 4 + 2A 2 + ie) 2 



(-2p°q° + 


P 


2 + 


q 


2 + s 4 + 2A 2 )(-2pV + 


IpI 


2 + 


q 


2 - si + S2S3/A&. + 2A 2 ) 


[(IP 




q |) 2 + 2A 2 ] [(|p| 


+ 


|q|) 2 - 


-2A 2 ] 



2p°g° + 2|p| 


|q 


~«4 (IP 




q|) 2 + 2A 2 


2p°g° - 2 


Pi 


|q 


-S4 (IP 


+ 


q|) 2 + 2A 2 



81 + s 4 - s 2 s 3 /Mj c , 
2|p||q| 

Now we will give some explanations about Eq. ( I7T1) . The results for k = q — p 
and = g + p + (rji — 772 )P have the same structures. We have changed the sign 
of the imaginary part of the pole for k = q + p + (771 — 77 2 )P by taking the variable 
transformation p° — > —p° . When k = q — p we will take X P ^(+P°> IpI) i n Eq. ()7ip 
and Sj (i = 1, . . . , 4) are defined by 

si = P 2 + g 2 + 4(77 1 2 -77 1 )P 2 + 2(2^ -l)(p-P + g-P), 

s 2 = p 2 - g 2 + 27/i(p • P - q ■ P) , 

s 3 = p 2 -g 2 + 2(77 1 -l)(p-P-g-P), 

s A = p 2 + q 2 -M^, 

while when k = q + p + (7/1 — 77 2 )P we will take x P ^K~ P°i IpI) i n Eq. flTTj) and 
(i = 1, . . . , 4) are defined by 

si = p 2 + q 2 -P 2 , 

S2 = P 2 ~q 2 - (2771 - 1)P 2 -2(77! " l)p-P-2r] iq -P, 

s 3 = p 2 -q 2 + (2r ll -l)P 2 ~2r ]l p-P-2(T ]l -l)q-P, 

s A = p 2 + g 2 -M^ + (2r7 1 -l) 2 P 2 - 2(277 1 -l)p-P + 2(277 1 -l)g-P. 

Now, we can substitute Eq. ( 1331) into Eq. (17B and complete the p°-integration by 
choosing proper contours. From the expression of f(p°) above we can see that the all 
the poles come from the BS wave function. This is because the denominator in f(p°) 
(neglecting the isospin violation, then 771 = 772) 

- 2p°q° + |p| 2 + |q| 2 + s A + 2A 2 = (p° - g ) 2 + 2A 2 - M 2 . 

is positive when A > M%* / V%, which is satisfied in our case (see the discussion in the 
next section). The remaining contour integration over p° is straightforward and the 
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result reads 
d 4 p 



(27T) 



F(|k|) 2 A^(k, M K .)i 2 (k + 2q 2 r(k - 2?!)" 

,o=- Vl E-E 1 + ^ 2 /(^ )Lo 



± 



k=q—p 



k=q+p 



)}x p (i \p) 



d\p\\p\ 2 [^f(p° 



P °= V1 E-Ei 



0> 



-TflE—E-2. 



p°=ri2E-E2 

(1) (|P|) 



(72) 



where fi = £ 2 (£ - S a - E 2 )/[(E 1 + £ 2 )(£ + E x - E 2 )}, & = E l {E + E l + E 2 )/[(E X + 
E 2 ){E + E x - E 2 )\ , £ 3 = E 2 {E + E x + E 2 )/[(Ei + E 2 ){E -E x + E 2 )\ , and £ 4 = 
Ei(E -Ei- E 2 )/[(Ei + E 2 )(E - E x + E 2 )] . If 771 = ^2 we have = ^4 and £ 2 = 
then for the iso- vector tttc final state, the decay width is zero. 

Once we have obtained the BS wave function of the ground state x p {p) (t ne numer- 
ical calculation will be carried out in the next section), we will take the wave function 
as input to calculate the decay amplitudes in Eqs. (|65l) and (|69|) . 



5 Numerical analysis and results 

The cutoff A in our model is not a free parameter in principle. It contains the informa- 
tion about the non-point interaction due to the structures of hadrons. In Ref . [7] , the 
cutoff for the interaction of KKp is taken to be rather large (about 3.18 GeV in our 
notation). On the other hand, in the study of baryons in the quark-diquark picture, 
the cutoff in the form factors associated with the diquark-gluon-diquark interaction is 
taken to be about 1.27 GeV [25]. In this work, we shall treat the cutoff A in the form 
factors as a parameter varying in a much wider range (0.8,4.8) GeV, in which we will 
try to search for possible solutions of the KK bound states. 

Let us first solve the BS bound state equation fl56|) numerically. We discretize the 
integral equation (I5U1) into a matrix eigenvalue equation by the Gaussian quadrature 
method. For each pair of trial values of the cutoff A and the binding energy E^ of the 
KK system (which is defined as E\> = E — mi — m 2 ), we will obtain all the eigenvalues 
of this eigenvalue equation. The eigenvalue closest to 1.0 for a pair of A and E^ 
will be selected out and called "the-trial-eigenvalue" . Fixing a value of the cutoff A 
and varying the binding energy E\> (from to —100 MeV) we will obtain a series of 
"the-trial-eigenvalue" s. For some (not all) values of the cutoff, we will find that the 
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Table 1: For the iso-scalar KK system, there are five regions of the cutoff A . In each 
region, for any value of the cutoff, the series of "the-trial-eigenvalue" s cross over the 
exact eigenvalue 1.0 (at certain binding energy E\> £ [—1, —99] MeV). 



E h ( MeV) 


A (GeV) 


-1 


1.1360 


2.0793 


2.7352 


3.5453 


4.7633 


-99 


1.2162 


2.0979 


2.7444 


3.5524 


4.7697 



Table 2: The decay widths (r^ ) corresponding to the five cutoff-regions when the 
binding energy E\> « —20 MeV. 



A (GeV) 


1.1700 


2.0862 


2.7385 


3.5479 


4.7656 


T (MeV) 


0.671 


4.758 


7.516 


11.102 


11.949 



corresponding series cross over i.o 13 in the range of E\> £ (0, —100) MeV. The task is 
then to find out all these cutoff values (which are, in fact, some continuous regions). 

In searching for the possible solutions in the iso-scalar channel of the KK system 
and its contribution to /o(980) (I G (J PC ) = + (0 ++ )), we find several regions of the 
cutoff. The results are listed in Table [U If we want to identify the iso-scalar bound 
state of the KK system with the observed /o(980) , we should let the binding energy 
approximately be —20 MeV. For the above five regions, this is equivalent to set the 
cutoff to be 

A « 1.1700, 2.0862, 2.7385, 3.5479, 4.7656 GeV, (73) 

respectively. From Eqs. ( 1651) and (l6"2l) . taking the corresponding BS wave function 
as input, the decay width of the iso-scalar KK to tttt can be obtained. The results 
are listed in Table [21 From the PDG's review [26] the full width of /o(980) is T = 
40 to 100 MeV (tttt dominant). Therefore, the results in Table [2] show that /o(980) can 
not be completely the KK iso-scalar bound state. 

Now, let us turn to the iso-vector KK bound state. We find the following two 
regions of the cutoff in this case: 



A £ (2.1160,2.2213) and (4.4998,4.5147) GeV. 
H That is to say, e.g., from 0.99 to 1.01 . 
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(74) 



The corresponding decay widths of the KK system into i\r\ with the binding energy 
E h = -20 MeV are given by 

r(. 7 = 1} = 1.329 MeV, A = 2.1590 GeV, (75) 
T^ 1 ) = 0.031 MeV, A = 4.5056 GeV. (76) 

The full width of a (980) is T = 50 to 100 MeV (nrj dominant) [26]. Although the KK 
iso-vector bound state does contribute to a (980) (I G (J PC ) = 1~(0 ++ )), just as in the 
case of / (980), a (980) can not be completely the KK iso-vector bound state. 

6 Conclusions and discussions 

In this paper we derive the BS equation for the KK system, study the possible bound 
states of this system, and calculate their decay widths in the BS formalism. In our 
model, we have used the ladder approximation. This approximation has been ques- 
tioned and is found not to be a good one in some models where higher order graphs give 
even more important contribution than the ladder graph [18-21]. However, in our case, 
we have shown explicitly that crossed-ladder graphs are suppressed greatly comparing 
with the ladder graphs due to the large masses of the exchanged particles. This makes 
the ladder approximation be legitimate in our model. In addition, based on the fact 
that the KK system is weakly bound, we have used the instantaneous approximation 
in the BS equation, in which the energy exchange between the constituent particles 
is neglected. Since the constituent particles and the exchanged particles in the KK 
system are not point-like, we introduce a form factor including a cutoff A which reflects 
the effects of the structure of these particles. Since A is controlled by non-perturbative 
QCD and can not be determined at present, we let it vary in a reasonable range within 
which we try to find possible bound states of the KK system. 

From the calculating results we find that there exist bound states of the KK system. 
Unfortunately, we can not determine the binding energy uniquely. The binding energy 
depends on the value of the cutoff A . For the iso-scalar KK system, we find five cutoff 
regions in which the solutions (with the binding energy E\> G (0, —100) MeV) to the 
ground state of the BS equation can be found (in unit of GeV): 

A ~ (1.136,1.216), (2.079,2.098), (2.735,2.744), (3.545,3.552), (4.763,4.770). 
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From these results, we can see that, except for the first interval, these regions are very 
narrow. For the iso- vector case, we find two regions (in unit of GeV), 

A ~ (2.1160, 2.2213) , (4.4998, 4.5147) . 

How to fix the cutoff (then the binding energy can be predicted), which is equivalent 
to how to determine the finite size effects of hadrons in the calculation, is beyond the 
scope of this paper. If we treat the binding energy as an input (E\> = —20 MeV), we 
find that the corresponding BS wave function gives too small decay widths, i.e. 

r^= 0) = 0.671, 4.758, 7.516, 11.102, 11.949 MeV, 

corresponding to the five cutoff regions, respectively. For the iso-vector case we have 

r^ =1) = 1.329, 0.0305 MeV, 

corresponding to the two cutoff regions, respectively. 

The authors in Ref. [10] concluded that the model with the one-meson-exchange 
potential from chiral dynamics, which is also used in this work, is sufficient to bind the 
KK system into a molecule which has the same mass and decay width as those of the 
iso-scalar / (980). From our calculation, however, we find that even these [KK) bound 
states could contribute to the observed scalar particles, the portion should be small 0. 
We prefer to draw the conclusion that there may be some more important structures 
besides the KK molecule in the observed overpopulated scalar particles (e.g. /o(980) 
and do (980)). Obviously, to resolve this problem further investigations are required. 
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